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We present a theory of the interplane conductivity of bilayer high temperature superconductors, 
focusing on the effect of quanta! and thermal fluctuations on the oscillator strengths of the superfluid 
stiffness and the bilayer plasmon. We find that the opening of the superconducting gap and estab- 
lishment of superconducting phase coherence each lead to redistribution of spectral weight over wide 
energy scales. The factor-of-two relation between the superfluid stiffness and the change below Tc in 
the oscillator strength of the absorptive part of the conductivity previously derived for single-layer 
systems, is found to be substantially modified in bilayer systems. PACS: 74.20-z,74.25.Gz,78.20.-e 



I. INTRODUCTION 

The interlayer ('c-axis') conductivity of high temper- 
ature superconductors is an important and long stand- 
ing problem. Experimental results iQ-^ have seemed to 
many workers |^-|^ to be sharply at variance with con- 
ventional understanding and to imply the existence of 
radically new physics. Other workers, conversely, have 
argued that many aspects of the results can be under- 
stood in a straightforward manner . Especially inter- 
esting have been apparent violations of the Ferrel-Glover- 
Tinkham sum rule |ll| relating the superfluid stiff'ness 
to changes in the absorptive part of the conductivity as 
temperature is reduced below the transition temperature 
Tc. 

Recently lofFe and one of us |^,^ have argued that 
the interlayer conductivity is a theoretically simple object 
(basically the convolution of two in-plane Green func- 
tions) and is therefore a sensitive probe of in-plane scat- 
tering rates and of the quantal and thermal phase fluctu- 
ations characterizing the superconducting state. A num- 
ber of predictions were made, some of which appear to 
agree with experiment and some of which do not |^ . The 
results reported in [T2[[r^ had a crucial limitation: the 
equations were derived for a 'single- layer' system such as 
La2-xSrxCuOi, whereas most (but not all) of the exper- 
imentally studied systems (including Bi2Sr2Ca2Cu20s 
and Y Ba2Cu307-s) have a 'bilayer' structure, with a 
unit cell containing two superconducting Cu02 planes 
coupled to each other more strongly than to the planes 
in adjoining unit cells. The new feature introduced by 
the bilayer structure is "local field corrections" : applica- 
tion of a uniform field can lead to a non-uniform charge 
distribution within a unit cell, which in turn causes inter- 
nal fields affecting the motion of charges. This leads to 
phenomena not found in single plane systems; for exam- 
ple, the bilayer plasmon feature observed and discussed 
by van der Marel and others jsj. Interest in this fea- 



ture was recently increased by the observation that 
the bilayer plasmon frequency may provide information 
about the in-plane electronic compressibility, a quantity 
of great theoretical interest not easily accessible by other 
techniques. 



In this paper we generalize the treatment of [|l2 13 
to the bilayer case. We provide a simple and physically 
transparent treatment of the c-axis conductivity in the 
limit (appropriate for the high temperature superconduc- 
tors) in which the interplane coupling is weak relative to 
in-plane energy scales. Our treatment includes phonon, 
bilayer plasmon and quasiparticle absorption. Our re- 
sults provide a justification for previously proposed phe- 
nomenological oscillator models and allow us to deter- 
mine the interplay between bilayer plasmon features and 
interlayer phase coherence. Our methods may easily be 
generalized to more complicated situations such as the 
three and four layer structures found in other high- Tc 
materials, but this generalization is not given here. 

The rest of this paper is organized as follows. In Sec II 
we present the formalism; Sec III gives results calculated 
in the absence of phonons; Sec IV discusses the spectral 
weight and sum-rule analysis. In Sec V we extend our 
treatment to incorporate phonons (relevant for some bi- 
layer materials) and finally in Sec VI we summarize our 
conclusions and discuss the applications to experiment. 



II. FORMALISM 

A. Fundamental Equations 

We study the bilayer system shown in Fig 1 in which 
each unit cell contains two conducting planes separated 
by a distance di and coupled by a hopping ti. The dis- 
tance between a plane in one unit cell and the closest 
conducting plane in another unit cell is ^2, so the lat- 
tice parameter in the interplane direction is d = di -I- c?2 • 
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Planes separated by a distance ^2 are coupled by a hop- 
ping ^2- We neglect further neighbor hoppings, although 
these can be easily added at the cost of increased com- 
plexity of our equations. In the high Tc context, t2 << ti 
(but our results are valid for any ratio ^2/^1) and both 
ti and t2 depend strongly on in-plane momentum, being 
maximal for momenta in the (0, tt) region of the zone and 
minimal for momenta near the zone diagonals (±7r, ±7r). 
We will usually not write the momentum dependence ex- 
plicitly. 



d, 




' applied 



Fig. 1 Geometry considered in present paper. Sliown is one 
unit cell (chosen as two close planes) and part of the next unit cell. 
The figure also shows the applied spatially uniform electric field 
{Eappiied), the charge build-up on planes (represented as + + ++ 

and and the resultant induced fields Ei„ji 2 which in turn 

affect the charge fiow. 

We refer to the two planes in one unit cell by the index 
I — a,b and to the region between two planes separated 
by di as region 1 and the region between two planes sep- 
arated by d2 as region 2; we label the unit cells by the 
index i. We take all planes to be identical and neglect 
all interplane couplings except for the hoppings and the 
internal electric fields induced by nonuniform charge dis- 
tributions. We allow for the possibility that the planes 
are at different electrochemical potential /i. The Hamil- 
tonian describing the system is then 

H ^ ^ ^171— plane (1) 
i,l 

where i7i„_p;ane (which we will not need to specify) de- 
scribes the electronic physics within a CuO^ plane. We 
shall study the properties of this Hamiltonian by a per- 
turbation expansion in t\ and i2- The dimensionless 
parameter is ti^2/ Ein^piane where Em-piane is the in- 
plane density of states or inverse of some other impor- 
tant in-plane local energy scale. This approach has been 
shown to agree with results obtained by other means in a 
number of contexts, including coupled Luttinger liquids 
| p5| , p^ and semiconductor heterostructures p7[ . 

We are interested in optical experiments^™ which 
may be thought of as involving the application to the sys- 
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tern of a weak spatially uniform transverse electric field 
of magnitude Et directed perpendicular to the planes. 
The experimentally determined quantity is the bilayer 
conductivity (JbUayer , which is the coefficient relating the 
applied electric field to the spatial average of the current. 
The applied field leads to an electrochemical potential 
fii^i on each plane which has three contributions: from 
the applied electric field, from fields generated by build- 
up of charge on particular planes (shown as Eindi,2 in Fig 
1), and from changes in the in-plane chemical potential 
due to changes in the in-plane density. We have 



(2) 



where Ri^i is the position vector of the plane in the in- 
terplane direction and the zeros of charge density n and 
of chemical potential /i have been defined to correspond 
to the states of the planes in equilibrium. eVind is the 
electrochemical potential due to electric fields produced 
by charge build up and = ed^/dn is inverse of the 
the exact in-plane density-density correlation function of 
Hin-piane- The factors of e arise from converting particle 
densities to charge densities. 

The spatially varying chemical potential leads to in- 
terplane electrical currents described by operators such 
as 



^'^ ^i,apa 



Ci,bpa- — H.C.) (3) 



and therefore to interplane charge build-up, for which 
we must solve self consistently. In the present simple 
situation, application of a uniform electric field leads to 
two independent densities Ua and rif, and two indepen- 
dent chemical potential differences, eVi = fXi^a — P^i,b and 
eV2 — /ij-1.6 — Mi. a- Combining the continuity equation 
for the current, the Maxwell equation and the density, 
and evaluating the currents to leading nontrivial order 
in ti,2 and Et leads to an expression for UbUayer- This 
expression is most conveniently written in terms of con- 
ductivities (Ji (CT2) appropriate to a 'single-layer' material 
consisting of an infinite stack of identical planes all sepa- 
rated by distance di ((^2) and coupled by hopping ti (12) 
and is 



^bilayer 



ai(72 - iuj{(Tidi + (J2d2)/C 
(Tid2 + o'2di — iuj/C 



(4) 



with di^2 — di.2/{di + (^2). Here the "blockade parame- 
ter" 
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e e'^did2 ' 



(5) 



where e is the 'background' dielectric function due to 
non-electronic degrees of freedom. In this paper we shall 
take e to be constant except in Sec V where we note in- 
clude the effects of phonons (important in the optical ab- 
sorption of some high- Tc materials) by using an e(w) with 
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appropriate frequency dependence in the expressions for 
o'biiayer- The constant C expresses the blocking effects 
arising because charge which flows onto a plane via the 
strong link (large conductivity) will take a long time to 
flow off via the weak link: if the driving frequency is low, 
charge buildup will therefore occur, inhibiting additional 
motion of charge across the strong link. 

Eq 1^ reproduces all of the obvious limits correctly: as 
^ it reduces to u'^aayer = (72 + di/ ai so uuiayer is 
dominated by smaller of the conductivities as expected, 
while ii u) — CIm ^<ti(w)c?2 + 'y2{'^)dij then a 'bilayer 
plasmon' pole occurs (damped, of course, by the dissipa- 
tive part of the conductivity). If di = d2 and cti = (T2 
then the system becomes effectively single-layered and 

Eq I shows auiayer = . 

The calculation of the constituent conductivities ai^2 is 
given in and relevant results will be recalled in the 

next subsection. We note here that for consistency they 
(and x^^) must be calculated to leading nontrivial order 
in the interplane hoppings ii,2. is therefore a single- 
plane quantity and (Ji^2 may therefore be expressed in 
terms of convolutions of two-dimensional in-plane Green 
functions. If higher order expressions are used then for 
example, exchange interaction contributions must be in- 
cluded in C and further changes to crtuayer will occur. 

We see that the frequency dependence of the bilayer 
conductivity is complicated and depends on the value 
of C and on the magnitudes and frequency dependences 
of the individual conductivities. In general <Ti,iiayer ex- 
hibits three regimes: a high frequency regime in which 
'^bilayer = cTidi + (T2'^2 ! 8- low frequency regime in which 
(^bilayer ~ o'2 / ^2 and a broad crossover regime with char- 
acteristic scale 



uj* ^ C cri{uj*)d2 + cr2{w*)di 



(6) 



which depends on the conductivities. If in the super- 
conducting state, Lu* < 2A then the scale oj* becomes 
identical to the bilayer plasmon frequency tObUayer and 
near WbUayer we have 



a{uj ^ LJbilayer < 2A) — 



'^Pbilayer 



^ — ^bilayer + iS 
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defining the strength npbuayer of the bilayer plasmon ab- 
sorption. 



B. Constituent conductivities 

The calculation of the constituent conductivities is dis- 
cussed at length in |l^,|l^. Here we briefly recall key 
results and needed formulae. The conductivities are 
given by correlation functions of current operators such 
as ji above and involve expectation values of the form 
t'^ip) < cf^pit)cj,p{t)c'^p{t')ci^p{t') > . To leading order in 



t, correlations between operators on different planes van- 
ish so the expression may be written as the sum of two 
terms, one involving < cf^p{t)ci^p{t') >< Cj'p(t')cj^p{t) > 
(i.e. the product of two 'normal' in-plane Green func- 
tions G{p,t — t')) and one involving < cfp{t)cfp{t') >< 
^3,p{^')^3,p{^) ^ (i-^- the product of two 'anomalous' in- 
plane Green functions F{p, t — t')). However, the anoma- 
lous Green function involves the superconducting order 
parameter which has a phase which we denote by </> . The 
product of anomalous Green functions on planes i and j 
therefore involves the factor e'^'^**^'"'*''"'^^^'' '* (times a 
short ranged function of r, t which depends on the de- 
tails of the interplane hopping and the underlying en- 
ergy scales of the superconductivity) and must be aver- 
aged over an ensemble describing the phase fluctuations. 
Rcfs ll|Jl| showed that in in the case of interest here, 
these effects may be accounted for by multiplying the 
F — contribution to cr by a Debye- Waller factor a 
which is unity for a mean-field BCS superconductor with 
no fluctuations, may be reduced from unity by quantal or 
thermal fluctuations about an ordered state, and which 
becomes very small if there is no long range phase order. 
We follow refs ||l^,|l^ in assuming that the pseudogap 
state is characterized by a conventional superconducting 
gap but no interplane phase coherence. 
Thus (z/ = 1, 2 labels planes) 
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with the diamagnetic contribution given by 

(-G{p, u^)G{p, J^) -f aF{p, Lu'jF{p, ^1)) (9) 
and the paramagnetic contribution given by 

= Ae'^d^Tj" I -^U{pf 

(g(p, Wn + wl)G(p, J J -t- aF{p, w„ + i^'n)F{p, ^1)) ■ (10) 

where G, F are the exact normal and anomalous Green 
functions corresponding to Hm-piane) ■ 
The w — > limit is cr — s- — with 

Ps,, = Sae'd,Tj2 I -^UipfFip,u;'jFip,u;'j (11) 

while the usual 'f-sum rule' arguments ||l^ yield (ct^^^ is 
the absorptive part of the conductivity) 



K ^ 



2dw 



TT ./0+ 



2dLJ 



a^-^\Lj) (12) 



(note that in the flrst equality the integral is only over 
one half of the delta function at a; = 0). 
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In the high- Tc materials the anisotropy of t{p) is such 
that the interplane conductivity is dominated by the 
"corner" regions around (0, tt) and so we foUow refs 
| p2| , p^ and neglect both the angular variation of the gap 
and of ti^2- In this approximation cti and a2 have the 
same frequency dependence and differ only by a prefac- 
tor involving the square of the relevant hopping. 

In high-Tc materials, the normal state c-axis conduc- 
tivity is characterized by a very broad Drude-like absorp- 
tion, corresponding to an in-plane Green function (in the 
"corner region" ) characterized by a very large, essentially 
frequency independent scattering rate. We therefore fol- 
low refs and use this 'dirty-limit' form to compute 
the conductivities. We note that there is a large and 
growing literature on changes to the Green function as 
the temperature is changed through the superconduct- 
ing transition [l9|-[2l|]. The implication of these changes 
for the c-axis conductivity has been studied by but 
because our main interest here is in the new features in- 
troduced by the bilayer structure and because there is no 
consensus on the physical origin or mathematical form of 
the superconductivity-induced changes, we do not con- 
sider them here. Further, we shall be interested mainly 
in three situations-the normal state, at a temperature 
well above the 'pseudogap formation temperature', the 
T — > limit in the superconducting state, and temper- 
atures well below the pseudogap scale and near to 
i.e. Tc < r << A. Thus we may neglect the tempera- 
ture, except as it influences the value of the Debye- Waller 
parameter a. We therefore have 



a«M = -(7o,.e(|L^|-2A) 



(13) 



duj 



^(2) ^ + ao,.sgn(^)e(2A - |c.|) (14) 

2\uj\ 

uj' juj' - \uj\) + a^A'^ dJ 



duj 



+fo, 1/3(771(^)0(1^1 — 2A) 

■'^1+^ uj'{lu' -\Lu\) + a^A^ 



a)|-A 



Vuj'^ -a^Ja^ - {uj'- 



^2 \UJ\ 



with (To,!/ the normal state (neither superconductivity nor 
any gap) conductivity, frequency independent because we 
have taken the dirty limit. We note that because ti 2 dif- 
fer, so also may the quantal fluctuation parameters ai,2- 
The considerations of |l^ suggest that a is dominated 
by in-plane fluctuations, so may not differ much between 
the two links, so in the rest of this paper we set ai = a2- 
The single-layer superfluid stiffnesses following from 
these expressions are 



where pf^^ is the superfluid stiffness following from the 
assumption of full phase coherence(aj/ = 1). 



III. CALCULATED CONDUCTIVITY 



In this section we evaluate the formulas derived in 
the previous section. The fundamental result was Eq 
^, which expressed the conductivity (JbUayer of a bilayer 
system in terms of the conductivities a 1,2 of effective 
'single-layer' systems corresponding to the two interplane 
spacings of the bilayer and a "blockade parameter" C ex- 
pressing interplane interaction effects. We use the normal 
state (A = 0) value of the characteristic frequency scale 
ijj* defined by Eq |, 



UJ 



Cao,id2 1 



(^ogdi 

(^Q,ld2 



(16) 



in our discussion henceforth. It is most convenient to 
express C in terms of this normal state value of oj* . 

The important dimensionless parameters are the ratio 
of normal state conductivities h = ctq, 2/fo,i < 1; i^* / A 
and the Debye-Waller factors 011.2 introduced above Eq 
^. To simplify the presentation of our results, we de- 
fine conductivity units such that cro,i — 1 and frequency 
units such that A = 1. For definiten^s we ch(X)se the 
normalized interplane distances to be di ~ 0.4, ^2 — 0.6 
so w* — 0.6cro,iC(l -I- 26/3). and set ai = a2- 
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Fig. 2 Real part of normal state conductivity scaled by 
0-0,1 {(^i + ^^2) for h = 1,0.5,0.25,0.075,0.01 (from top to bot- 
tom). The bilayer frequency oj* = 11.48,9.91,8.8,8.39 for h = 
0.5, 0.25, 0.075, 0.01, respectively. The inset shows imaginary part 
of the conductivity for b = 0.075. 

Fig 2 shows the real part of the normal state con- 
ductivity for the five values b = 1,0.5,0.25,0.075,0.01. 
The suppression of the low frequency conductivity by the 
blockade effect is evident, as is the gradual crossover to 
the high frequency isolated layers value. The curves have 
been scaled by ctq, i(di -I- 6(^2) so that all have the same 
high frequency limit. The inset shows the imaginary part 
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for b = 0.075. In the crossover regime, the blockade effect 
is seen to lead to out of phase response. 

We now consider the superconductivity induced 
changes. The top panel of Fig 3 shows the result of 
evaluating Eq |^ in the 'single-layer' (b = 1) fully phase 
coherent {a„ — 1) case. The opening of the supercon- 
ducting gap suppresses the real part of the conductivity 
for frequencies below 2A and changes the form somewhat 
for oj ^ 2 A. The establishing of superconducting phase 
coherence leads to a divergent low frequency response 
characterized by the superfluid stiffness Ps,biiayer- The 
oscillator strength in the superfluid response is shown in 
the top panel as a shaded rectangle. The f-sum rule ar- 
guments discussed at length in the next section imply 
that in the fully phase coherent {a^, = 1) case, the area 
lost in the absorptive part of a due to the opening of the 
superconducting gap is transferred to the superfluid re- 
sponse. It is apparent from the figure that the area in the 
shaded rectangle is approximately equal to the 'missing' 
area and we have verified numerically that the areas are 
equal: J^^ duj[a{A = 0) - cr(A)] = 7rps,biZayer/2. 



a= 1 



uj < 2A region by the opening of the gap will go into the 
bilayer plasmon absorption instead of into the superfluid 
delta function. 
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Fig. 3 Real part of conductivity for 'single-layer' (b = 1) system. 
Top panel corresponds to fully phase coherent (a^ = 1) while the 
bottom panel to =0.1 superconducting state. The dashed line 
in both corresponds to the normal state. The area of the shaded 
rectangle equals nps/2 in each case. 

If 6 < 1 then the situation is much more involved. In 
particular because the form of the conductivity in the 
regime w ~ w* depends sensitively on the interplay be- 
tween C and (Ti.2, the superconductivity induced changes 
in cri.2 will lead to large changes in (TbUayer- In addi- 
tion, some of the oscillator strength eliminated from the 



Fig. 4 Real part of conductivity for bilayer (6 = 0.075) sys- 
tem for 1^;* = 0.4A. Top panel corresponds to fully phase coherent 
(ai, = 1) and bottom panel to = 0.1 superconducting state. 
The dashed line in both corresponds to the normal state conduc- 
tivity. The area of the two shaded rectangles in each panel equals 
nps/2 and irpbUayfir foi' rectangles centered at a; = and (^bUayer^ 
respectively. 



The top panels of Figs 4, 5 and 6 show the normal 
and superconducting state conductivity in the fully phase 
coherent (a^ = 1) limit for three representative values 
of u!*: 0.4A, 3A and lOA, respectively in the strongly 
anisotropic limit b = 0.075. The oscillator strengths 
in the superfluid response and (if it is inside the gap) 
the bilayer plasmon are shown as shaded rectangles. For 
UJ* = 0.4A, the bilayer plasmon lies within the supercon- 
ducting gap. The resulting absorption is a delta function 
at the marked frequency, with an intensity correspond- 
ing to an integrated area equal to that of the rectangular 
box shown. For uj* — 3A the bilayer plasmon lies just 
above the gap, visible as a sharp feature at the gap edge 
and the remainder of the superconducting conductivity 
is slightly suppressed over a wide frequency range. For 
UJ* — lOA the bilayer plasmon feature is evident only 
as a very broad absorption at frequencies that are much 
larger than shown here. 
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Fig. 5 Real part of conductivity for bilayer (6 = 0.075) sys- 
tem for u)* = 3A. Top panel corresponds to fully phase coherent 
{civ = 1) and bottom panel to oti, = 0.1 superconducting state. 
The dashed line in both corresponds to the normal state conduc- 
tivity. The area of the shaded rectangle at oj = in each panel 
equals irps/s while the area of the rectangle at uibUayer in the sec- 
ond panel equals iTpi,iiayer- 



The bottom panels of Figs 4, 5 and 6 show the ef- 
fect of phase fLuctuations, displaying the superconducting 
curves for the same lo* and h values as the respective top 
panels, but with a,^ = 0.1. For lo > 2 A the difference be- 
tween the normal and the superconducting conductivity 
increases as phase fluctuations become more important. 
Both the strength and the frequency of the bilayer plas- 
mon feature depend strongly on the value of the fluctua- 
tion parameter. For lu* = 3A, IDA, the bilayer plasmon 
moves below the gap as seen in Figs 5 and 6. For the 
pseudogap case {a^, — 0) both the bilayer plasmon and 
the superconducting delta function are absent and the 
Lu > 2 A conductivity is roughly the same as the = 0.1 
case. We can also study the conductivity for different 
values of ai and a2 and it is worth noting that in the 
case when a2 = and ai ^ 0, it is possible to get a bi- 
layer plasmon feature though the superconducting delta 
function is absent. 



CO* = 1 A 



0.4 



- 


1 ' 1 ' 

a=1 

- 








a = 0.1 


1 


1 



2 4 

CO/ A 



Fig. 6 Real part of conductivity for bilayer {b = 0.075) sys- 
tem for 1^;* = lOA. Top panel corresponds to fully phase coherent 
(ai, = 1) and bottom panel to = 0.1 superconducting state. 
The dashed line in both corresponds to the normal state conduc- 
tivity. The area of the shaded rectangle at = in each panel 
equals nps/s while the area of the rectangle at uituayer in the sec- 
ond panel equals TTptiiayer- 

Fig 7 plot the bilayer plasmon frequency oJuiayen the 
spectral weight in the bila^yer plcLsmon PbUayer 

(defined 

by Eq 0) and the spectral weight in the superconducting 
delta function(at a; = 0) as a function of ai = ai 
for 6 = 0.1 and uj* = 0.4 A. The two spectral weights 
vary linearly in the fluctuation parameter for the given 
value of w* and hence their ratio is independent of the 
ay value. 
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Fig. 7 The bilayer plasmon frequency uihUaycr^ the spectral 
weight in the bilayer plasmon PhUayer (defined by Eq |^ and the 
spectral weight in the superconducting delta function (at tj = 0) 
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ps as a function of a\ = 02 for h = 0.075 and ui* = 0.4A. 



IV. SUPERFLUID STIFFNESS AND SUM RULES 

At very low frequencies in the superconducting state 
one has 



^Ps,bilayer 



yer 



(17) 



Inspection of Eq Q shows that 

Ps,lPs,2 



Ps^bilayer 



Ps,ld2 + Ps,2dl 



= TTCTo 1 A- 



baia2 



aid2 + ba2di 



(18) 



where the first equation apphes to all bilayer systems and 
the second follows from the specific assumptions we have 
made. One question of current experimental interest is 
the relation between the superfluid stiffness and changes 
in conductivity as the temperature is reduced (a) below 
the 'pseudogap' scale at which the gap opens and (b) be- 
low Tc at which phase coherence is established. In Refs 
these relations were established for the 'single- 
layer' case. The numerical results presented in the pre- 
vious section show that differences occur in the bilayer 
case. To analyze this issue more precisely, we note that 
the dirty-limit model analyzed in this paper should be 
viewed as arising from a model with a very large but fi- 
nite scattering rate T in the limit {Q, ^bUayen A) << F. 
The standard sum rule derivations are based on analysis 
of the f2/r — > 00 limit. However, one may obtain sum 
rules for the superconductivity and pseudogap induced 
changes in a without considering this limit. We define 
the change in the spectral weight as A is increased from 
A = by 



SK{n,A,a)^ / — [(t(A, a) - (7(A 
Jo 



0)]. (19) 



This quantity remains finite in the limit 17 — > 00, rj/F << 
1. The values we obtain with our F — > 00 limit are accu- 
rate up to terms of relative order {^bUayer, A) /F. 

It is also useful to consider the change in the spectral 
weight excluding the superfluid response : we define 

SK+{n, A,a)^ / — [a(A,a)-CT(A = 0)] 
Jo+ 

= SK{fl,A,a) - ps{A,a) (20) 

and the ratio of change in SK^ with a to the superfluid 
stiffness given by 



R{n) = 



_ SK+{Q,A,a) - SK+{Q, A,0) 



P. (A, a) 



(21) 



Now Refs [|2|,|3j showed that for a 'single-layer' (s-1) 
system, the change in the spectral weight as defined by 
Eq[l|is 

dKs-i[U = 00, A, a) = , (22) 

where pf^^lA) is defined by Eq |l^. Thus the change 
with a in total spectral weight at fixed A is 



SKs-iin ^00, A, a) -dK,-iin^ 00, A,0) 
apf^^A) 



(23) 



and the change with a in the > spectral weight as 
defined by Eq pOl is 



SK+^,_i{n = cx), A,a) - 5K+^s-ii^ = oo, A,0) 



ap 



BCS 



(A) pr'(A,«) 



(24) 



Use of Eq. 21 gives the value of Rs-i{^ — 00) — —1/2. In 
other words, if the superconducting gap appears without 
phase coherence, the oscillator strength decreases by an 
amount related to the 'BCS' superfluid stiffness, essen- 
tially because the conductivity in the region less than the 
gap is suppressed and no additional oscillator strength 
appears in the superfluid response. If phase coherence 
is now turned on, the total oscillator strength and the 
superfluid stiffness increase, while the spectral weight in 
the > conductivity decreases. Comparison of Eq 14 
and Eq 23 shows that in the 'single-layer' case the ratio 
between these changes is 2:1. 

Applying these arguments to the bilayer case shows 
that 

1 



SKbiiayeri^ = OO, A, a^,) = -- 2_^(1 - ai,)7rao,i,d^A 



(25) 



The simple factor-of-two relation between the phase- 
coherence-induced change in SK^ and the superfluid stiff- 
ness does not occur in the bilayer system essentially be- 
cause when the phase coherence parameter is varied, the 
strengths of both the bilayer plasmon feature and the su- 
perfluid stiffness vary. To see what the relation is, we plot 
in Fig 8 the ratio RbUayeri^ = 00) as a function of bi- 
layer anisotropy b {RbUayeri^ — 00) being independent of 
ai = a2 = a). We see that the ratio increases monoton- 
ically as b is decreased from the single-layer value b = 1, 
and changes sign at a di/d2 dependent value of b ^ 0.2. 
Thus unlike in the single-layer case, where Ps increased by 
twice the decrease in SK^{Q = 00), in the bilayer case 
the increase is generically greater than 26K+{n — 00) 
and indeed for extreme anisotropy both SKj^(J} — 00) 
and Ps increase as a is increased from zero. 
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II 

a 



DC 




Fig. 8 Ratio of the change in 5K+ on onset of phase coherence 
to the superfiuid stiffness defined by Eq RbUayeri^ = oo) 
a function of the anisotropy parameter 6. Note that the ratio is 
independent of the value of a and that the 6 = 1 value corresponds 
to the single-layer case. 

Fig 9 plots Rbiiayer as E fuiiction of the cut-off il. The 
— oo value is indicated by an arrow. We have verified 
that the calculated quantity docs indeed converge to the 
correct D, —^ oo value, but as can be seen, the convergence 
is very slow. 



CO* = 0.4 A 



a 



DC 
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0.5 





2 3 
Q./ A 



Fig. 9 Ratio of the change in on onset of phase coher- 

ence to the superfiuid stiffness defined by Eq RbUayer plotted 
as a function of cut-off Q for a bilayer {b = 0.075) system with 
uj* = 0.4A. Solid line: fully phase coherent (a^ = 1) supercon- 
ducting state; dotted line: «^ =0.1 superconducting state. Arrow: 

Fig 10 shows the changes in spectral weight for the case 
of a large difference in interplane hopping (6 = 0.075) and 
for Lu* — 0.4A as a function of the cut-off frequency fl. 
For a = 1 (top panel), we expect conservation of the to- 
tal spectral weight while for a = 0.1 (bottom panel) we 
expect the D, oo value to be non-zero as given by Eq 
p5| . A remarkably slow convergence of the change in the 
spectral weight to its 17 ^ oo value is evident. We ver- 
ify the f2 — > oo values numerically by plotting in Fig 11, 
SKhiiayer as a function of inverse cut-off frequency. We 
have verified that the l/fi ^ limit matches the value 
given by Eq 25. Considerable caution must be exercised 
in the experimental investigations of changes in the spec- 



tral weight because small differences persisting over wide 
frequency ranges may lead to appreciable contributions 
to sum rules. 



0.4 A 




Fig. 10 SKbiiayer (defined by Eq as a function of cut-off 
frequency f2 for bilayer (fe = 0.075) system for u)* = 0.4A. Solid 
line corresponds to fully phase coherent (a^ = 1) and dotted line 
to olu — 0.1 superconducting state. 

Q)" = 0.4 A 




0.02 0.03 

A/n 

Fig. 11 SKbiiayer '^^ a function of inverse cut-off frequency Q 
demonstrating convergence to correct sum-rule value. The 1/Q = 
value correctly gives the f2 — > oo value as given by Eq Top panel 
corresponds to fully phase coherent {a^ = 1) and bottom panel to 
ai, =0.1 superconducting state. 
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V. INCLUSION OF PHONONS 



Due to the proximity of the bilayer feature to the opti- 
cal phonons in high- Tc materials like Bi2Sr2Ca2Cu20% 
and Y Ba2Cu^OT-Si it is of interest to study the inter- 
play of phonons and the bilayer plasmon. In the first 
part of this paper we had taken £i = £2 = £00 = const. 
To incorporate phonon modes we include in our analysis 
frequency dependent dielectric functions for each layer: 
Si{lo) and £2(w). We obtain for interplane bilayer con- 
ductivity, 



phonon 
bilayer 



(to) 



- |^(o-i(£2 - d2) + (J2{ei - di)) 



aid2 + 0-2^1 



4tt 



[eid2 + £2^1) 



16 



7^(£l£2 - (£1^2 + £2^1)) 



(71^2 + <J2dl - |j(£lrf2 + £2^1 



(26) 
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This more general formula reduces to Eq ^ with a value 
of C = in/e on choosing £1 = £2 = constant (inclusion 
of the compressibility term x leads to more complicated 
formulae). 

To study the qualitative effect of including phonons we 
consider the simplest possible case of €2 = 1 and 



£l(w) 



(£0 ^ £00)1.^^^ 



(27) 



where LUp is the frequency of the phonon mode and 7 the 
broadening. 

Fig. 12 displays the effects of adding, to the situation 
(full phase coherence (a = 1) and bilayer plasmon in- 
side the gap) shown in the top panel of Fig 4, a phonon 
with a frequency greater than the bilayer plasmon fre- 
quency. For orientation, the top panel shows the bilayer 
plasmon part of the electronic absorption in the absence 
of phonons (calculated from Eq.^ with £1 = £2 = 1) , 
and the phonon absorption in the absence of electrons 
(calculated from Eq. ^with cri = 0-2 = 0). The bilayer 
plasmon was represented in Fig 4 as a rectangle and is 
shown here with a Lorentzian broadening. The phonon 
feature is rendered optically active and shifted up from 
the phonon frequency ujp by bilayer effects. The lower 
panel of Fig 12 shows the full conductivity. The elec- 
tronic continuum contribution at a; > 2A is present but 
very difficult to perceive on the scale of this plot. It is ev- 
ident that the coupling between the modes leads as usual 
to level repulsion, and that further, almost all of the os- 
cillator strength goes into the upper mode. The lower 
mode (shown in expanded view in the inset) is almost 
invisible. 



Fig. 12 Top panel plots the conductivity CTphononi^) witli 
ojp = 0.91 A obtained from Eq^by putting tlie electron conductiv- 
ities (71^2 = (dotted line) and the broadened bilayer conductivity 
o'bUayeri^^) with uj* = 0.4A and b = 0.075 (dot-dashed Une). The 
bottom panel plots o'jj';n^°"('*^) (Eq^). The inset shows the lower 
peak not visible on the scale of the plot. 
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Fig. 13 Top panel plots the conductivity dphononi'-^) with 
uip ~ 0.23 A obtained from Eq^by putting the electron conductiv- 
ities (Ti 2 = (dotted line) and the broadened bilayer conductivity 
utiiayeri'^) with UJ* = 0.4A and b = 0.075 (dot-dashed line). The 
bottom panel plots o'jj'in^"" ('^) (Eq p^). The inset shows the lower 
peak not visible on the scale of the plot. 
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Fig. 13 shows that roughly the same situation is ob- 
tained if the phonon starts out at a lower frequency than 
the bilayer plasmon. In this case the upper mode shifts 
by rather less, but the qualitative features are the same. 

We now briefly outline the effects of increasing phase 
fluctuations (decreasing a from unity) i.e increasing tem- 
perature. As can be seen from Figs 4-6, increasing 
phase fluctuations decreases the frequency and oscillator 
strength of the bilayer plasmon mode. Thus if the 'bare' 
phonon frequency is greater than the bilayer plasmon fre- 
quency (as in Fig 12), then relatively minor changes occur 
in the absorption spectrum as a is decreased. Essentially, 
the almost invisible lower absorption moves to the left 
and becomes a bit sharper which in turn results in the 
slight decrease in frequency and intensity of the upper 
absorption. 
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a = 0.1 
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CO/ A 
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Fig. 14 The figure plots <,';°;;°"(<.^) (Eq for a = 0.8 (top 
panel), a = 0.3 (panel two) and a = 0.1 (bottom panel) for the 
same set of parameters as in Fig 13 [uip ~ 0.23A lu* = 0.4A). 



On the other hand, if the a = 1 bilayer plasmon is 
at a higher frequency than the phonon, more drastic 
changes will occur as shown in Fig 14. Decreasing a 
slightly results in a decrease in frequency and intensity 
of the stronger upper absorption as shown in the top 
panel of Fig 14 for a = 0.8. As a is further decreased 
a the Ubiiayer peak in Fig 13 goes on moving to the left 
and eventually overlaps with the Uphonon peak at which 
point the lower absorption is almost invisible and the 
conductivity is as shown in the second panel of Fig 14 
for a = 0.3. On further decreasing a the (Tuiayer peak 
crosses over to the left of (Jphonon and the the conductiv- 
ity is as given by the bottom panel of Fig 14 for a = 0.1. 
Also note that the decrease in intensity of the upper peak 
is much stronger with the decrease in a as compared 
to the case when the 'bare' phonon frequency is greater 
than the bilayer plasmon frequency. From experimen- 
tal point of view, the behavior of the observed peaks as 
a function of temperature might allow us to distinguish 
between the two situations where the 'bare' phonon fre- 
quency is greater than or smaller than the bilayer plas- 
mon frequency. 

Finally, we note that the spectral weight in the super- 
fluid stiffness is unaffected by the phonons. 



CONCLUSIONS AND APPLICATIONS TO 
EXPERIMENT 



VI. 



We have extended the theory of refs ||12|,|13[ to the bi- 
layer situation of relevance to many experimentally stud- 
ied high-Tc materials. The crucial physics is local-field or 
"blockade" effects: the difference in hopping amplitudes 
characteristic of a bilayer structure leads to charge imbal- 
ances inside the unit cell; the electrochemical potentials 
due to these charge imbalances act to suppress the low 
frequency response to a uniform electric field, and lead 
to bilayer plasmon features in the absorption. For the 
physics of high-T the crucial question is the observabil- 
ity of the effects of thermal and quantal fluctuations of 
the phase of the superconducting order parameter, here 
parameterized by a 'Debye- Waller parameter a 

We find a low frequency suppression of the conductiv- 
ity in the normal (neither superconducting or pseudogap) 
state discussed in section III and shown in Fig 2. However 
this distinct signature is not apparent in the experimen- 
tal plots of the normal state conductivity. There are two 
possibilities: cither the scale is very high or the effect is 
masked by phonons. 

A second, generally valid qualitative result is that the 
simple factor of two relation between ps and the change 
with onset of phase coherence in the uj > oscillator 
strength which was found for single-layer systems, no 
longer applies for bilayer systems (cf Eq ^ and Fig 8); 
the change in ps is generically greater. This qualitative 
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behavior has been observed. 

A third important result is that the convergence of 
sum-rule integrals with frequency can be very slow; so 
caution should be exercised in applying sum-rule argu- 
ments to data. 

Further we showed that the coupling of phonons, phase 
fluctuations and the bilayer plasmon leads to complicated 
effects on the spectrum, which depend sensitively on pa- 
rameters suggesting that an unambiguous extraction of 
the bilayer plasmon frequency and spectral weight may 
be difficult. This is unfortunate, as the these quantities 
in principle carry information about the physically cru- 
cial phase fluctuation properties encoded in the Debye- 
Waller parameter a. We suggest however how the change 
is conductivity with increasing phase fluctuations (or in- 
creasing temperature) is expected to be different based 
on whether the phonon frequency is above or below the 
bilayer plasmon frequency. 
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